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We outline the construction of metastable de Sitter vacua of type IIB string theory. Our starting 
point is highly warped IIB compactifications with nontrivial NS and RR three-form fluxes. By 
incorporating known corrections to the superpotential from Euclidean D-brane instantons or gaugino 
condensation, one can make models with all moduli fixed, yielding a supersymmetric AdS vacuum. 
Inclusion of a small number of D3 branes in the resulting warped geometry allows one to uplift the 
AdS minimum and make it a metastable de Sitter ground state. The lifetime of our metastable de 
Sitter vacua is much greater than the cosmological timescale of 10 10 years. We also prove, under 
certain conditions, that the lifetime of dS space in string theory will always be shorter than the 
recurrence time. 
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I. INTRODUCTION 

There has recently been a great deal of interest in find- 
ing de Sitter (dS) vacua of supergravity and string the- 
ory. This is motivated in part by the desire to construct 
possible models for late-time cosmology (since a small 
positive cosmological constant seems to be required by 
recent data Q), and in part by more conceptual wor- 
ries that arise in the study of dS quantum gravity (as 
discussed, for instance, in 0). The no-go theorem of 
Q guarantees that such solutions cannot be obtained in 
string or M-theory by using only the lowest order terms in 
the lOd or lid supergravity action, but one expects that 
corrections to the leading order Lagrangian in the g s or 
a' expansion or inclusion of extended sources (branes) 
should improve the situation. Indeed, a careful discus- 
sion of how such additional sources (which are present in 
string theory) invalidate the no-go theorem for warped 
backgrounds and allow one to find highly warped com- 
pactifications appears in p|. Additional sources which 
violate the assumptions of the theorem were shown to 
yield dS vacua in noncritical string theory in JjJ. 

Here, we use our knowledge of quantum corrections 
and extended objects in string theory to argue that there 
are dS solutions of ordinary critical string theory. Our 
basic strategy is to first freeze all the moduli present 
in the compactification, while preserving supersymmetry. 
We then add extra effects that break supersymmetry in 
a controlled way and lift the minimum of the potential 
to a positive value, yielding dS space. To illustrate the 
construction we work in the specific context of IIB string 
theory compactificd on a Calabi-Yau (CY) manifold in 
the presence of flux. As described in Q such construc- 
tions allow one to fix the complex structure moduli, but 
not the Kahler moduli of the compactification. In par- 



ticular, to leading order in a' and g s , the Lagrangian 
possesses a no-scale structure which does not fix the over- 
all volume (we shall assume that this is the only Kahler 
modulus in the rest of this note; it is of course possible 
to construct explicit models which have this property). 
In order to achieve the first step of fixing all moduli, 
we therefore need to consider corrections which violate 
the no-scale structure. Here we focus on quantum non- 
perturbative corrections to the superpotential which are 
calculable and show that these can lead to supersymme- 
try preserving AdS vacua in which the volume modulus 
is fixed in a controlled manner. 

Having frozen all moduli we then introduce supersym- 
metry breaking by adding a few D3 branes in the com- 
pactification; this is a compactified version of the situa- 
tion discussed in 6] . The addition of D3 branes does not 
introduce additional moduli: their worldvolume scalars 
are frozen by a potential generated by the background 
fluxes |{J. Inclusion of anti-D3 branes in the absence 
of other corrections yields a run-away to infinite volume 
of the compact space (since the energy density in the 
D3 branes generates a tadpole for the volume modulus). 
However, we show that in the presence of the quantum 
corrections we have described, in a sufficiently warped 
background the D3 tension can be a small enough cor- 
rection to lift the formerly AdS vacuum to positive cos- 
mological constant, without destabilizing the minimum. 

The extent of supersymmetry breaking and also the re- 
sulting cosmological constant of the dS minimum can be 
varied in our construction, within a range, in two ways. 
One may vary the number of D3 branes which are in- 
troduced in the above manner, and one can also vary 
the warping in the compactification (by tuning the num- 
ber of flux quanta through various cycles). It is impor- 
tant to note that this corresponds to a freedom to tune 
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discrete parameters, so while fine-tuning is possible, one 
should not expect to be able to tune to arbitrarily high 
precision. Since there is still a vacuum at very large ra- 
dius with approximately zero energy (this is the Dine- 
Seiberg runaway vacuum 0), any dS minimum is only 
a false vacuum; but it is only destabilized by tunneling 
effects, and we argue that the lifetimes one can achieve 
are extremely long. In addition, we argue that under 
the assumption that the potential between the dS mini- 
mum and the Dine-Seiberg vacuum at infinity is positive 
(which will be true in any simple examples, since a single 
dS maximum is the only intervening critical point), the 
lifetime of the dS minimum is always shorter than the 
timescale for Poincare recurrences discussed in Q . 

While our emphasis in this paper is on dS vacua it 
is worth remarking that the first step of our construc- 
tion, which freezes all moduli, is of interest in its own 
right (for another recent approach to this problem, see 
(tj; other recent ideas about fixing the volume modulus, 
while leaving other moduli unfixed, appear in |lflj and 
references therein). In fact moduli stabilization has been 
an important open question in string theory and related 
phenomenology. We show that this can be achieved by 
putting together a few different effects, all of which are 
quite well understood by now. Preserving susy allows us 
to carry out this analysis with control. Once all moduli 
are frozen, susy breaking effects other than the introduc- 
tion of D3 branes alone can also be considered. Some 
of these, like the D3/D7 inflationary models of ^]|, or 
models with both branes and anti- branes [12j , are also of 
possible cosmological interest. 

A more complete discussion of the possible cosmologi- 
cal toy models that one can construct using combinations 
of the ingredients described in this paper, in the spirit of 
0, will appear in 

We should note that there has been a great deal of 
interesting work on constructing dS solutions in super- 
gravity and string theory. dS minima of 4d gauged super- 
gravities which do not as yet have a known string theory 
embedding appeared in [l5L Il6| , while dS compactifica- 
tions of gauged 6d supergravity appear in |17| . The work 
of [f| constructs dS vacua in supercritical string theory 
using many of the same ingredients which arise here. Fi- 
nally, the importance of using fluxes in the cosmological 
context was stressed in ^8|, where however the problem 
of moduli stabilization was left as a black box (related 
ideas appeared in [19j). Cosmology of the simplest flux 
compactifications, on the T 6 jZi orientifold pol l2ll |22| 
(whose gauged supergravity description is worked out in 
123), was recently investigated in |24| . 



II. FLUX COMPACTIFICATIONS OF IIB 
STRING THEORY, INCLUDING CORRECTIONS 

In this section, we briefly describe the required knowl- 
edge of flux compactifications of type IIB string theory. 
In section IIA we describe the models of 0, and in 



section IIB we enumerate various quantum corrections 
which can modify the superpotential and Kahler poten- 
tial used in 0. In section IIC, we show that incorporat- 
ing the generic corrections can yield (supersymmctric) 
AdS minima with all moduli stabilized. 



A. Calabi-Yau orientifolds with flux 

We start with F-theory [2{| compactified on an elliptic 
CY fourfold X. The F-theory fourfold is a useful way of 
encoding the data of a solution of type IIB string theory; 
the base manifold M of the fibration encodes the IIB 
geometry, while the variation of the complex structure r 
of the elliptic fiber describes the profile of the IIB axio- 
dilaton. In such a model, one has a tadpole condition 

Here T3 is the tension of a D3 brane, Nm is the net 
number of (D3 — D3) branes one has inserted filling the 
noncompact dimensions, and H3, F3 are the three- form 
fluxes in the IIB theory which arise in the NS and RR 
sector, respectively. As shown by Sen [2(|, in the absence 
of flux, it is always possible to deform such an F-theory 
model to a locus in moduli space where it can be thought 
of as an orientifold of a IIB Calabi-Yau compactification. 
For this reason, we will use the language of IIB orien- 
tifolds, with M being the Calabi-Yau threefold which is 
orientifolded. In this language, the term x ffl counts the 
negative D3-brane charge coming from the 03 planes and 
the induced D3 charge on D7 branes, while the terms on 
the right-hand side count the net D3 charge from trans- 
verse branes and fluxes in the CY manifold. As in Q , we 
will assume we are working with a model having only one 
Kahler modulus, so h}' x {M) = 1 (/i M (A) = 2, and one 
modulus is frozen in taking the F-theory limit, where one 
shrinks the elliptic fiber). Such models can be explicitly 
constructed, by e.g. using the examples of CY fourfolds 
in |27| or by explicitly constructing orientifolds of known 
CY threefolds with h 1 ' 1 = 1. 

In the presence of the nonzero fluxes, one generates a 
superpotential for the Calabi-Yau moduli, which follows 
from |28| (see also [2{J[3(J) and is of the form 

W = f G3ACI (2) 
Jm 

where G3 = F$ — tH^, with r the IIB axiodilaton. Com- 
bining this with the tree-level Kahler potential 

K = -31n[-i(p-p)]-ln[-i(r-f)]-ln[-i / OaH] 

Jm 

(3) 

where p is the single volume modulus (p = bj \/2+ie 4u ^; 
our conventions are as in £|), and using the standard 
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N = 1 supergravity formula for the potential, one finds 



V = e K \^g al D a WD b W -Z\W\ 

\ a,b 

e K (^f j DiW~DjW) (4) 



Here, a,b runs over all moduli fields, while i,j runs over 
all moduli fields except p; and we see that because p 
does not appear in (J2J, it cancels out of the potential 
energy (@J, leaving the positive semi-definite potential 
characteristic of no-scale models |3l| . 

One should use this potential as follows. Fix an in- 
tegral choice of #3,-^3 in H 3 (M, Z); then, the potential 
@ fixes the moduli at values where the resulting G3 
is imaginary self-dual (ISD). Supersymmetric solutions 
furthermore require G3 to be type (2,1) (more generally, 
G3 would have a (0,3) piece). Thus in supersymmetric 
solutions W = on the vacuum, while in the nonsuper- 
symmetric solutions, W = Wq, a constant which is de- 
termined by the (0,3) piece of G3. In generic solutions, 
the complex structure moduli of the F-theory fourfold (in 
IIB language, the complex structure moduli, the dilaton, 
and the moduli of D7 branes) arc completely fixed, leav- 
ing only the volume modulus p. The scale of the masses 
?7i for the moduli which are fixed is 
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(5) 



where R is the radius of the manifold (Im p scales like 
i? 4 ). In this approximation, R is unfixed. By tuning flux 
quanta, it is possible (at least in some cases) to fix g s at 
small values, though not arbitrarily small. 

There is one last point we will need to use in sec- 
tion III. The fluxes (and any transverse branes) serve 
as sources for a warp factor. Therefore, such models 
with branes and flux are generically warped compacti- 
fications. In fact, as shown in Q, following earlier work 
of |23, H2, IH, 0, |3f|, it is possible to construct mod- 
els with exponentially large warping. One can write the 
Einstein frame metric of the compactification as 

ds\ Q = e^rj^dx^dx" + e~ 2A ^~g mn (y)dy m dy n (6) 

with y coordinatizing the compact dimensions, and g mn 
the unwarped metric on M (so in the orientifold limit, 
it is a Calabi-Yau metric). Then it is shown in H (by 
compactifying the Klebanov-Strassler solution |3||) that 
one can construct models parametrized by flux integers 
M, K such that 



exp[-(27r^)/(3 5s M)] 



(7) 



with e being of order one at generic points. This means 
in particular that with reasonably small flux quanta, one 
can generate exponentially large ratios of scales in such 
models. 



In the following, we will assume that g s and the com- 
plex structure moduli have been fixed at the scale © 
by a suitable choice of flux, and concentrate on an ef- 
fective field theory for the volume modulus p. This is 
self-consistent, in that the final mass for p will be small 
compared to (J5J. One last comment. Light states could 
also arise from modes living in the throat region which 
experiences large warping. We assume here that any such 
excitations are gapped, as in [3(j. Then typically the p 
modulus, which has a Planck scale suppressed mass, will 
be much lighter than these excitations and we can neglect 
them as well in the low energy theory. 



B. Corrections to the no-scale models 

Here, we write down two known sources of corrections 
to the no-scale models, both parametrize possible correc- 
tions to the superpotential J2J- Then, in section IIC, we 
show that including either correction to the superpoten- 
tial yields supersymmetric models with AdS vacua. 

1) Witten has argued that in type IIB compactifica- 
tions of this type, there can be corrections to the super- 
potential coming from Euclidean D3 branes [37||. This 
happens when the fourfold X used for F-theory compact- 
ification admits divisors of arithmetic genus one, which 
project to four-cycles in the base M. In the presence of 
such instantons, there is a correction to the superpoten- 
tial which at large volume yields a new term 



Winst = T(zi) exp(27rip) 



(8) 



where T{zi) is a complex structure dependent one-loop 
determinant, and the leading exponential dependence 
comes from the action of a Euclidean D3 brane wrap- 
ping a four-cycle in M. Since the z\ and the dilaton are 
fixed by the fluxes at a scale ©, we can integrate them 
out and view JSj as simply providing a superpotential for 
the volume modulus. 

2) In general models of this sort, one finds (at special 
loci in the complex structure moduli space) non-Abelian 
gauge groups arising from geometric singularities in A, or 
in type IIB language, from stacks of D7 branes wrapping 
4-cycles in M. Assume that the fluxes have fixed one at 
a point in moduli space where this phenomenon occurs 
(examples appear in |38j , for instance) . Consider a stack 
of N c coincident branes. The 4d gauge coupling of the 
SU(N C ) Yang-Mills theory on such wrapped branes (we 
ignore the decoupled £7(1) factor) satisfies 



8tt 2 
9ym 



R A 

2tt — = 2?r Im p 

9s 



(9) 



Since the complex structure moduli of X are completely 
fixed, the D7 brane moduli (at least in cases where the 
4-cycle being wrapped has vanishing h , which are easy 
to arrange) are also fixed. Therefore, any charged mat- 
ter fields (which would create a Higgs branch for the D7 
gauge theory) have also been given a mass at a high-scale; 
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and the low-energy theory is pure N = 1 supersymmet- 
ric SU(N C ) gauge theory. This theory undergoes gluino 
condensation, which results in a nonperturbative super- 
potential 

Wgauge = A^ c = Ae^ (10) 

where Ajy c is the dynamical scale of the gauge theory, 
and the coefficient A is determined by the energy scale 
below which the the SQCD theory is valid (There are also 
threshold corrections in general, these contribute sub- 
leading effects.) We see that this leads to an exponential 
superpotential for p similar to the one above (but with a 
fractional multiple of p in the exponent, since the gaug- 
ino condensate looks like a fractional instanton effect in 
W). 

So effects 1) and 2) have rather similar consequences 
for our analysis; we will simply assume that there is 
an exponential superpotential for p at large volume. In 
our companion paper [l4j|. we investigate some interest- 
ing possibilities for cosmology if there are multiple non- 
Abelian gauge factors. Using the fourfolds in [23, it is 
easy to construct examples (with h l,1 {X) — 2) which 
could yield gauge groups of total rank up to ~ 30. The 
results of [3!j suggest that much larger ranks should be 
possible. 

One important comment is in order before we proceed. 
Besides corrections to the superpotential of the kind dis- 
cussed above, there are also corrections to the Kahler 
potential (see e.g. f° r a calculation of some lead- 
ing corrections). In our analysis we will ensure that the 
volume modulus is stabilized at values which are para- 
metrically large compared to the string scale. This makes 
our neglect of Kahler corrections self consistent. 



solved by turning on only flux, i.e. with no additional 
D3 branes. 

At a supersymmetric vacuum D p W — 0. We simplify 
things by setting the axion in the p modulus to zero, and 
letting p — ia. In addition we take A, a and Wo to be all 
real and Wq negative. The minimum then lies at 

DW = -> W = -Ae- a " er (l + \a(j cr ) (13) 

The potential, V = e K (G p P D p WI)pW - 3\W\ 2 ), at 
the minimum is negative and equal to 

VAds = (~3e K W 2 ) Ads = ~-^r (14) 

D (Tcr 

We see that we have stabilized the volume modulus while 
preserving supersymmetry. It is important to note that 
the AdS minimum is quite generic. Any corrections to 
the Kahler potential will still result in a susy minimum 
which solves ijTT?)) . 

A few comments are in order before we proceed. A 
controlled calculation requires that a ^> 1, this ensures 
that the supergravity approximation is valid and the cJ 
corrections to the Kahler potential are under control. It 
also requires that a a > 1 so that the contribution to 
the superpotential from a single (fractional) instanton is 
reliable. Generically, if the fluxes break supersymmetry, 
Wq ~ O(l), and these conditions will not be met. How- 
ever it is reasonable to expect that by tuning fluxes one 
can arrange so that Wo <C 1. In these circumstances we 
see from (|13|l that a a > 1. Taking a < 1, one can then 
ensure that a ^> 1, as required. 

As an illustrative example we consider Wo = — 10~ 4 , 
A = 1, a = 0.1. This results in a minimum at cr cr ~ 113. 



C. Supersymmetric AdS Vacua 

Here, we show that the corrections to the superpoten- 
tial considered above can stabilize the volume modulus, 
leading to a susy preserving AdS minimum. We perform 
an analysis of the vacuum structure just keeping the tree- 
level Kahler potential 

K = -3hx[-i(p-7)] (11) 

and a superpotential 

W = W + Ae lap . (12) 

Wo is a tree level contribution which arises from the 
fluxes. The exponential term arises from either of the 
two sources above, and the coefficient a can be deter- 
mined accordingly. In keeping with the fact that the 
complex structure moduli and the dilaton have received a 
mass (0 , we have set them equal to their VEVs and con- 
sider only the low-energy theory of the volume modulus. 
To avoid the need to worry about additional open-string 
moduli, we assume the tadpole condition has been 
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FIG. 1: Potential (multiplied by 10 15 ) for the case of expo- 
nential superpotential with Wo = — 1CP 4 , A = 1, o = 0.1. 
There is an AdS minimum. 

Another possibility to get a minimum at large vol- 
ume is to consider a situation where the fluxes preserve 
susy, and the superpotential involves multiple exponen- 
tial terms, i.e. "racetrack potentials" for the stabilization 
of p |4l| . Such a superpotential could arise from multiple 
stacks of seven branes wrapping four cycles which cannot 
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be deformed into each other in a susy preserving manner 
[14) . In this case by tuning the ranks of the gauge groups 
appropriately one can obtain a parametrically large value 
of a at the minimum. 

Finally it is plausible that in some cases, there could 
be some other moduli which have not yet been fixed in 
this approximation (if, for instance, the fluxes were non- 
generic and left some of the complex structure moduli of 
the Calabi-Yau manifold unfixed). Since the dilaton is 
already fixed by fluxes and the radial modulus is fixed 
due to non-perturbative exponential terms in the super- 
potential, one might find interesting and realistic cosmol- 
ogy driven by these other scalar fields from the compact- 
ified string theory. Here, however, we will neglect this 
possibility, and focus instead on another mechanism for 
uplifting the AdS vacuum to a dS vacuum. 



III. CONSTRUCTING DS VACUA 

In this section, we uplift the supersymmetric AdS 
vacua of section IIC to yield dS vacua of string theory. In 
section IIIA, we describe the new ingredient: D3 branes 
transverse to M. In section IIIB, we show that for rea- 
sonable choices of parameters, the inclusion of D3 branes 
yields dS vacua. 



D3 branes in ISD fluxes 



(|15|l arises because the ambient five-form flux adds a re- 
pulsive energy equal to the tension for an D3 brane, see 
[42j (above equation (3.9)). Note also that we are con- 
sidering solutions which meet the ISD condition, even in 
the presence of the additional flux required to insert the 
D3 brane. There are corrections to l|15|l but these are 
quadratic in the number of D3 that are added and are 
small. 

The important point is that due to the warping the 
addition of the D3 brane breaks supersymmetry by a 
very small amount. In general terms, we get a term in 
the potential which goes like 



SV = 



8D 



(Im p) 2 



(16) 



(the factor of 8 is added for later convenience). The coef- 
ficient D depends on the number of D3 branes and on the 
warp factor at the end of the throat. These parameters 
can be altered by discretely changing the total flux, and 
the fluxes which enter in J7J), respectively. This allows us 
to vary the coefficient D and the susy breaking in the sys- 
tem, while still keeping them small (More properly, since 
the flux can only be discretely tuned, D can be varied 
but not with arbitrary precision). We will see that by 
tuning the choice of D one can perturb the AdS vacua of 
IIC to produce dS vacua with a tunable cosmological con- 
stant. The vacua will clearly only be metastable, since 
all of the sources of energy we have introduced vanish as 
Im p — > 00 . 



In the tadpole condition Q), there is a contribution 
from both localized D3 branes and from fluxes. To find 
AdS vacua with no moduli in the previous section, we 
assumed that the condition was saturated by turning on 
fluxes in the compact manifold. Now, we assume that 
in fact we turn on too much flux, so that Q can only 
be satisfied by introducing one D3 brane. In the flux 
background determining a solution of the sort described 
in section II, the D3 does not have translational moduli; 
they are fixed by the ISD fluxes, which generate a poten- 
tial for the world volume scalars. This kind of situation 
was studied in |6j. 

Now, the tadpole is cancelled, but there is an extra bit 
of energy density from the "extra" flux and D3 brane. 
In fact, as in equation (73) of 6], one finds that the D3 
brane adds an additional energy: 



SV 



, a T 3 



1 



9s (Imp) ; 



(15) 



with ao the warp factor at the location of the D3 brane. 
As described in |6j, in the presence of ISD fluxes of the 
sort characterizing the Klebanov-Strassler throat, any 
anti-D3 branes are driven to the end of the throat, where 
the warp factor is minimized. It follows from (jTJ) that 
the value of ao is exponentially small, and hence in suit- 
able models the inclusion of a D3 adds to the potential 
an exponentially suppressed term. The prefactor of 2 in 



B. Uplifting AdS vacua to dS vacua 

We now add to the potential a term of the form D/a 3 , 
as explained above. For suitable choices of D, the AdS 
minimum will become a dS minimum, but the rest of 
the potential does not change too much. There is one 
new important feature, however: there is a dS maximum 
separating the dS minimum from the vanishing potential 
at infinity. The potential is: 



V 



aAe 



2a 2 



aaAe- aa + W + Ae 



D 

7^3 



By fine-tuning D, it is easy to have the dS minimum very 
close to zero. For the model Wq = — 1CP 4 , A = 1, a = 0.1 
D = 3x 1(T 9 we find the potential (multiplied by 10 15 ): 
Note, if one does not require the minimum to be so 
close to zero, D does not have to be fine-tuned so pre- 
cisely. A dS minimum is obtained as long as D lies within 
a range, eventually disappearing for large enough D. If 
one does fine tune to get the minimum very close to 
zero, the resulting potentials are quite steep around the 
dS minimum. In this circumstance, the new term ba- 
sically uplifts the potential without changing the shape 
too much around the minimum, so the p field acquires a 
surprisingly large mass (relative to the final value of the 
cosmological constant). 
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FIG. 2: Potential (multiplied by 10 15 ) for the case of ex- 
ponential superpotential and including a correction with 
D = 3 x 10~ 9 which uplifts the AdS minimum to a dS mini- 
mum. 

It is important to mention that the value of the volume 
modulus shifts only slightly in going from the AdS mini- 
mum to the new dS minimum. This means if the volume 
was large in the AdS minimum to begin with, it will con- 
tinue to be large in the new dS minimum, guaranteeing 
that our approximations are valid. 

If one wants to use this potential to describe the 
present stage of acceleration of the universe, one needs 
to fine-tune the value of the potential in dS minimum to 
be Vo ~ 10~ 120 in units of Planck density. In principle, 
one could achieve it, e.g., by fine tuning D. However, 
the tuning we can really do by varying the fluxes etc. in 
the microscopic string theory is limited, though it may 
be possible to tune quite well if there are enough three- 
cycles in M. 



IV. HOW STABLE IS THE DS VACUUM? 

The radial modulus a = Im p has a kinetic term 
-^2 {da) 2 which follows from the Kahler potential 
For cosmological purposes it is convenient to switch to 

the canonical variable tp = J~^\n.a = ^/|ln(Imp), 

which has a kinetic term ^(dtp) 2 . In what follows we 
will use the field tp and it should not be confused with 
the dilaton, 4>. 



A. General theory 

The dS vacuum state ipo corresponding to the local 
minimum of the potential with Vo > is metastable. 
Therefore it may decay, and then the universe will roll to- 
wards large values of the field ip and decompactify. Here 
we would like to address two important questions: 

1) Do our dS vacua survive for a large number of 
Planck times? For instance, if we fine tune to get a small 
cosmological constant, is the dS vacuum sufficiently sta- 



ble to survive during the 10 10 years of the cosmologi- 
cal evolution? If the answer is positive, one can use the 
dS minimum for the phenomenological description of the 
current stage of acceleration (late-time inflation) of the 
universe. 

2) Is the typical decay time of the dS vacuum longer 
or shorter than the recurrence time t r ~ e s °, where 
S = is the dS entropy 0]? If the decay time 

is longer than t r ~ e s °, one may need to address the 
issues about the consistency of the stringy description of 
dS space raised in 0, IE @ • 

We will argue that the lifetime of the dS vacuum in our 
models is not too short and not too long: it is extremely 
large in Planck times (in particular, one can easily make 
models which live longer than the cosmological timescale 
~ 10 10 years), and it is much shorter than the recurrence 
time t r ~ e s °. 

In order to analyse this issue we will remember, fol- 
lowing Coleman and De Luccia |44| , basic features of the 
tunneling theory taking into account gravitational effects. 

To describe tunneling from a local minimum at tp = ipo 
one should consider an 0(4)-invariant Euclidean space- 
time with the metric 

ds 2 = ch 2 + b 2 {T){dijj 2 + sin 2 VdOf) . (17) 

The scalar field tp and the Euclidean scale factor (three- 
sphere radius) b(r) obey the equations of motion 

p" + 3 b -<p> = V v , h " = -^ + V ), (18) 

where primes denote derivatives with respect to r. (We 
use the system of units M p = 1.) 

These equations have several instanton solutions 
(tp(r), b(r)). The simplest of them are the 0(5) invari- 
ant four-spheres one obtains when the field tp sits at one 
of the extrema of its potential, and b(r) = H^ 1 sin Ht. 
Here H 2 = y, and V(tp) corresponds to one of the ex- 
trema. In our case, there are two trivial solutions of this 
type. One of them describes time-independent field cor- 
responding to the minimum of the effective potential at 
tp = tpo, with Vo = V(tpo). Another one is related to the 
maximum of the potential at tp = tpi, with V± = V(tpi). 

Coleman-De Luccia (CDL) instantons are more com- 
plicated. They describe the field tp(r) beginning in a 
vicinity of the false vacuum tpo at t = 0, and reaching 
some constant value tpf > tpx at r = r/, where b(rf) = 0. 
It is tempting to interpret CDL instantons as the tunnel- 
ing trajectories interpolating between the different vacua 
of the theory. However, one should be careful with this 
interpretation because the trajectories tp{r) for CDL in- 
stantons do not begin exactly in the metastable minimum 
tpo and do not end exactly in the absolute minimum of 
the effective potential. We will discuss this issue later. 

According to [4^| , the tunneling probability is given by 

P(cp) = e - s ^ +So , (19) 
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where S((p) is the Euclidean action for the tunneling tra- 
jectory (<ys(r), &(t)), and So = S(ip ) is the Euclidean 
action for the initial configuration ip — ipo- Note that 
S(ip) in Eq. (|19|l for the tunneling probability is the in- 
tegral over the whole instanton solution, rather that the 
integral over its half providing the tunneling amplitude. 
The tunneling action is given by 



%>) = J d* Xy /g(-±R- 



1 



(dip) 2 + V(<p) 



(20) 



In d = 4 the trace of the Einstein equation is R = (dip) 2 + 
4V(<p). Therefore the total action can be represented by 
an integral of V(<p): 

S(<p) = -f d i x^V(^) = -27r 2 J f drb 3 (r)V( V (T)) . 

(21) 

The Euclidean action calculated for the false vacuum 
dS solution ip = ipo is given by 



(22) 



Similarly, for the dS maximum ip — ip\ one has S\ = 

24tt 2 

This action for dS space So has a simple sign-reversal 
relation to the entropy of de Sitter space Sq: 



So — —Sq 



2Att 2 
~~Vo 



(23) 



Therefore the decay time of the metastable dS vacuum 
^decay ~ P^ 1 ' W) can be represented in the following way: 



+ . _ P S( V )+S _ f S( V ) 



(24) 



The semiclassical approximation is applicable only for 
| *S* (v 3 ) I ^ 1- Eq. (|21(l implies that for the tunneling 
through the barrier with V(p) > (which is the case 
for the tunneling from dS space to Minkowski space in 
our model) the action S(tp) is always negative, S(jp) < 0. 
This means that the decay time of dS space to Minkowski 
space is exponentially smaller than the recurrence time 
t r . The existence of the runaway vacuum at oo in field 
space with zero energy is a standard feature of all string 
theories Q ■ We conclude that the problems related to the 
decay time exceeding the recurrence time t r 0, HI 0] will 
not appear in the simplest string theory models, where a 
single dS maximum will separate the dS minimum from 
oo. 

Now that we found the upper bound on the tunneling 
time, we will try to estimate the tunneling time in our 
model using some particular instanton solutions. In gen- 
eral, it is very difficult to find analytical solutions for the 
CDL instantons and calculate the tunneling probability. 
We will investigate this problem using two different ap- 
proaches, which are, in a certain sense, opposite to each 
other: the thin-wall approximation and the no-wall ap- 
proximation. 



B. Thin- wall approximation 

Let us assume that one can describe the present stage 
of the acceleration of the universe by our model, so that 
Vq ~ 10 -120 in Planck units (more generally, the analysis 
of this section will work for Vq very small compared to 
the starting AdS cosmological constant, which can be ar- 
ranged by tuning) . This is a hundred orders of magnitude 
smaller than the height of the barrier V\ in our model. In 
Minkowski space, the condition Vq <C V\ usually means 
that the thin-wall approximation is applicable. Let us 
check whether one can use it in our case. 

In application to our scenario (tunneling from dS space 
with vacuum energy Vo > to Minkowski space) one can 
represent the results of |44| in the following useful form: 



P — exp — 



S(yo) 



(1 + (4V /3T2))s 



(25) 



Here T is not temperature but tension of the bubble wall; 
in our case 



T = 



d>p ^2V(p) 



(26) 



Eq. H25|) confirms our general conclusion that the sup- 
pression of the tunneling is always smaller than e -8 ^ -*, 
and the decay time is shorter than t r . 

There are two limiting cases of special interest: Vo 3> 
T 2 and Vo -C T 2 . The meaning of these two condi- 
tions becomes clear if we restore M p in these inequalities: 
V M 2 > T 2 and V M 2 <C T 2 . If we turn off the gravita- 
tional interactions (M p — > oo), one has VqM 2 ^> T 2 , and 
we obtain the well known result for the tunneling in the 
thin- wall approximation in Minkowski space: 



P : 



exp 



27tt 2 T 4 



(27) 



From Eq. I|26ll one finds that T ~ A^y^Tj where Aip 
is the typical width of dS maximum. This means that 
Vo 3> T 2 (in units M p — 1), and one can ignore the 
gravitational effects in the thin-wall approximation only 
if 



(28) 



One can easily check that for our model this condition 
is not satisfied, so we must take gravitational effects into 
account and study an opposite limit Vo -C T 2 . Then in 
the first approximation one simply has 



P w exp (-S(^o)) = exp 



24tt' 



exp(-10 122 ) 



(29) 

Thus, for all practical purposes our dS vacuum is com- 
pletely stable. 
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On the other hand, if one expands Eq. i|25|) in powers 
of V /T 2 , one finds 



P « exp (— S(ip a )) x exp 



647T" 
p2 



(30) 



For the sub-Planckian tension T <C 1 one finds that the 
probability of the tunneling is exponentially larger than 
P « exp (— S(yjo))- The decay time of dS space due to 
CDL instantons i^ecay ^ s mucn smaller than the recur- 



rence time t r 
general result: 



exp (— S(ipo)), in agreement with our 



CDL 
decay 



t r exp 



6471^ 
p2 



(31) 



In the thin- wall approximation (for Vq -C T 2 ) the ra- 



n (fo 
0, 



dius of the bubble is Rb — 4/T |44|, whereas the thick- 
ness of the wall is approximately ARb = T/V%. This 
means that the thin-wall approximation taking gravity 
into account is valid if ARb <C Rb, he. T 2 <C Vi- 



C. "No- wall" approximation: Hawking-Moss 
instanton and stochastic approach to tunneling 

Now we will return to the simplest instanton solutions, 
sitting at the dS minimum ip = tpo and the dS maximum 
ip = (p\. According to Hawking and Moss (HM) [45|. they 
describe tunneling through the barrier, with the proba- 
bility of tunneling suppressed as 



p = e -S{ V i)+S{ V0 )) 



I • x : i 



24tt 2 



24^ 



(32) 



This result may seem rather controversial because the in- 
stanton solution <p(r) = tfx does not interpolate between 
the stable vacuum and the false vacuum. 

In fact, as we already mentioned, the last problem ap- 
pears for all CDL solutions as well. These solutions never 
begin exactly in the false vacuum, so how could they be 
considered interpolating solutions? This problem disap- 
pears in the thin- wall limit, but it shows up very clearly 
in numerical calculations going beyond the thin wall ap- 
proximation. 

There are several different ways in which one can ad- 
dress this issue. In it was noticed that instead of 
considering an exactly constant solution <p = ipx one 
may consider a configuration that coincides with this 
solution everywhere except a small vicinity of the end- 
point at Tf = ttH^ 1 . The action involves integra- 
tion J dTb 3 (T)L(g tiU ,ip). Since the scale factor b(r) = 
H^ 1 sin Ht vanishes at Tf — ttH^ 1 , one can make very 
strong modifications of the solution ip = ipx ha a small 
vicinity of Tf without changing the action in a signifi- 
cant way. If the action changes by less than O(l), then 
each such configuration can be used for the description of 
tunneling despite the fact that it is not a true solution of 
classical equations of motion near r/. This allows one to 



bend the solution (p = ipx (or the CDL solution) in such 
a way that it begins to interpolate between the different 
vacua. 

One may also patch HM instantons and CDL instan- 
tons to each other in several different ways, which may 
provide their alternative interpretation and justify the 
results obtained by using CDL and HM instantons, see 
e -g- E3- 1-3 OH- More recently a possible interpreta- 
tion of CDL and HM instantons was suggested in [50| . 
Significant progress was achieved by Gen and Sasaki in 
the development of a consistent Hamiltonian approach to 
tunneling with gravity |5l|. 

But these methods do not address another problem 
with the Hawking-Moss tunneling: Since the instanton 
(p = (px is exactly homogeneous, it seems to describe a 
homogeneous tunneling in the whole universe, which is 
impossible. One can circumvent this problem by claiming 
that the whole universe is reduced to the interior of a 
single causal patch of size H . However, this approach, 
which is often used for the description of eternal dS space 
(which does not decay), may be less useful in applications 
to inflationary cosmology since it loses information about 
the whole universe except for a small part of size H . 

The most intuitively transparent description of the 
(ne arly ) homogeneous tunneling was provided in [H^, IH, 
|5J, log in the context of the stochastic approach to infla- 
tion. 

One may consider quantum fluctuations of a light 
scalar field tp with m 2 = V" <C H 2 = V/3. During 
each time interval St = H^ 1 this scalar field experiences 
quantum jumps with the wavelength ~ H^ 1 and with 
a typical amplitude Sip = H/2ir. Then the wavelength 
of these fluctuations grows exponentially. As a result, 
quantum fluctuations lead to a local change of the am- 
plitude of the field ip which looks homogeneous on the 
horizon scale if -1 . From the point of view of a local ob- 
server, this process looks like a Brownian motion of the 
homogeneous scalar field. If the potential has a dS mini- 
mum at ipo with m <C H, then eventually the probability 
distribution to find the field with the value tp becomes 
time-independent: |53. 153 . l54l l55 | 



( \ 

P(tp, ip ) - exp ^"^y J ■ ex P 



This result was obtained without any considerations 
based on the Euclidean approach to quantum gravity. 
It provides a simple interpretation of the Hawking-Moss 
tunneling. During inflation, long wavelength perturba- 
tions of the scalar field freeze on top of each other and 
form complicated configurations, which, however, look 
almost homogeneous on the horizon scale H^ 1 . If origi- 
nally the whole universe was in a state (po , the scalar field 
starts wandering around, and eventually it reaches the 
local maximum of the effective potential at ip = <px- Ac- 
cording to ll-i-il) , the probability of this event is suppressed 

by exp ^— ^p? — h 2 ^-^j. As soon as the field ip reaches 

the top of the effective potential, it may fall down to 



9 



another minimum, because it looks nearly homogeneous 
on a scale of horizon, and gradients of the field ip are 
not strong enough to pull it back to ip . The probability 
of this process is given by the Hawking-Moss expression 
(|32|l . However, this is not a homogeneous tunneling all 
over the universe, but rather a Brownian motion, which 
looks homogeneous on the scale H~ x [53|. 

One can take a logarithm of the probability distribu- 
tion l|33(l and find that the entropy of the state with the 
field stochastically moving in the potential V(ip) is given 
by 



24^ 2 



(34) 



This is the simplest derivation of dS entropy that does 
not involve ambiguous Euclidean calculations. It is valid 
even for the states outside the dS minimum, as long as 
the condition \m 2 \ <C H 2 (i.e. \V"\ <C V) remains valid. 

This result allows one to obtain a simple interpreta- 
tion of the HM tunneling, proposed in [I]]]. Indeed, Eq. 
(|32|l has the standard thermodynamic form describing 
the probability of thermal fluctuations 



(35) 



A similar thermodynamic approach was recently devel- 
oped in a series of papers by Susskind et al § . 

Since the entropy S(<pi) — 24n 2 /Vi is positive, one 
finds another confirmation of our general result: 



t. 



HM _ -Si+So 

decay e 



t r exp 



24tt- 

' Vx 



< t r 



(36) 



condition coincides with the inflationary slow-roll condi- 
tion \V"\< V. 

Thus we are coming to the conclusion that the HM 
tunneling, and the thermodynamic approach discussed 
above, are most efficient for the description of the tun- 
neling in the inflationary universe, where their validity 
has been firmly established by the stochastic approach 
to inflation. On the other hand, in the situations where 
the potentials are very thin, A<f> -C 1, one should use 
the CDL instantons, and the thin- wall approximation is 
valid. 

This result has a very simple interpretation. CDL in- 
stantons describe tunneling through the barrier. Mean- 
while, HM instantons in the inflationary regime can be 
interpreted in terms of the Brownian motion, when the 
field slowly climbs to the top of the barrier due to ac- 
cumulation of quantum fluctuations with the Hawking 
temperature H/2tt. If the barrier is very wide, it is eas- 
ier to climb the barrier rather that to tunnel through it. 
In this case the HM tunneling prevails and the stochas- 
tic/thermodynamic description of this process in terms of 
dS entropy is very useful. If the barrier is very thin, it is 
easier to tunnel, and CDL instantons are more efficient. 

For the simple models with the parameters given in 
our paper one has 3T 2 < 8Vi, and the tunneling occurs 
mainly due to CDL instantons. Even in this case the 
stochastic/thermodynamic approach may remain useful: 
If this approach is valid not only in the context of infla- 
tionary cosmology with \V"\ <C V but also in the situa- 
tions with \V"\ > V, it provides a simple upper bound 
on the decay time of dS space, idecay < e 



S0po)-S0pi) 



On the other hand, in our case S\ <C So, so in the first 
approximation one finds, as before, that i^ecay ~ t r ~ 

1Q 122 

e 

Note that the quasi-homogeneous HM tunneling and 
the CDL tunneling correspond to two different processes; 
depending on the potential, one of these processes may 
happen much faster than another. Let us compare the 
rate of the quasi-homogeneous HM tunneling with the 
rate of tunneling due to CDL instantons in the thin wall 
approximation for T <C 1 : 



jHM 

decay 
+CDL 

decay 



exp 



8 - { T*~Vx 



(37) 



This shows that t™ ay 
nates) for 3T 2 > 8Vi 



< t 



CDL 

decay 



(HM tunneling domi- 



Meanwhile, in the opposite case, 
3T 2 < 8Vi, the tunneling occurs due to CDL instantons. 
This is consistent with our estimate of validity of the thin 
wall approximation, T 2 <C V±. 

It is useful to represent this result in a different way. 
Let us use the estimate T ~ A(p\/Vi, where Aip is the 
typical width of dS maximum. This estimate implies that 
3T 2 > 8Vi and t™ ay < t^ y if Atp > 1, i.e. if the width 
of the maximum is much greater than the Planck mass. 
For the potentials with the nearly Minkowski minima this 



V. DISCUSSION 

It has been a difficult problem to construct realistic 
cosmologies from string theory, as long as the moduli 
fields are not frozen. In this note, we have seen that it is 
possible to stabilize all moduli in a controlled manner in 
the general setting of compactifications with flux. This 
opens up a promising arena for the construction of string 
cosmologies. 

More specifically, we have seen that it should be pos- 
sible to construct metastable dS vacua in the general 
framework of Q, by including anti-branes |(| and incor- 
porating non-perturbative corrections to the superpoten- 
tial from D3 instantons |37j or low-energy gauge dynam- 
ics. For cosmology, it might be more interesting to in- 
clude a more complicated sector at the end of the warped 
throat. For instance, the D3/D7 inflationary models of 
[IT| . or models with both branes and anti-branes [l2| . 
could be of interest in this context. It should also be 
possible to make similar constructions in the heterotic 
theory, perhaps using the models of 56] as a starting 
point (and then incorporating a mechanism to stabilize 
the dilaton in a controlled manner). 

It is worth summarizing why we believe the construc- 
tions in this paper are reliable. Our analysis was car- 
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ried out in the framework of supergravity. There are two 
kinds of corrections to this analysis, in the g s and a' ex- 
pansions. By tuning the fluxes we have argued that both 
of these can be controlled. For example, it is reasonable 
to believe that g s can be stabilized at a value of order 
g s ~ 0.1 and the volume a 3> 1 (for the volume this was 
discussed in section II C). Note that for reliability we do 
not require that the dilaton and volume be made arbi- 
trarily small or big; in fact this is not possible since the 
fluxes can only be tuned discretely. We only require that 
these moduli take appropriately small or big values, and 
this can be achieved, especially if M has enough three- 
cycles, yielding many possible choices of flux background. 

We were also able to prove, in section IV, that our 
dS minima are short lived compared to the timescale for 
Poincare recurrences. This generalizes to any construc- 
tion where the dS minimum is separated from the Dine- 
Seiberg run-away vacuum (which is ubiquitous in string 
theory) by a potential which remains non-negative; the 
simplest controlled examples in string theory will have 
this feature. One can also imagine more complicated 
shapes of the potential between the dS minimum and in- 
finity, which include some intervening AdS critical points. 
It is natural to wonder if a more general statement can 
be obtained that would apply in these cases as well. 



Finally, it would be interesting to discover string the- 
ory models which naturally incorporate both > 60 e- 
foldings of early universe inflation, and a late-time cos- 
mology in agreement with the most recent data Q. (To 
show that it is possible to obtain a small enough value 
of the dS cosmological constant for late-time cosmology, 
one would have to demonstrate that an idea along the 
lines of |18| can be implemented in this context.) Some 
further steps towards making more realistic cosmological 
toy-models inspired by string theory, in the same general 
framework as this note, will be presented in [l4j . 
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